QUADRATIC HARNESSES FROM GENERALIZED BETA INTEGRALS 



WLODEK BRYC 



Abstract. We use generalized beta integrals to construct examples of Markov processes with 
linear regressions, and quadratic second conditional moments. 



1. Introduction 



O 

(N 

1.1. Quadratic harnesses. In [BM W07] the authors consider square-integrable stochastic pro- 
cesses on (0, oo) such that for all t, s > 0, 

<N '• (1.1) E(X t ) = 0, E(X t X s ) = min{t, s}, 

E(X t |J r SiU ) is a linear function of X S ,X U , and Var[Xi | J 7 ,^] is a quadratic function of X S ,X U . 
P_! Here, J 7 ,^ is the two-sided a-field generated by {X r : r G (0, s] U [u, oo)}. Then for all s < t < u, 

en \ (II. ip implies that 

> . 

C3 ■ / x s it — t t — s 

1.2 E (X t J" = X s + X u , 

u — s u — s 

which is sometimes referred to as a harness condition, see [MY05] . While there are numerous 
■ examples of harnesses that include all integrable Levy processes ( [JP881 (2.8)]), the assumption 
of quadratic conditional variance is more restrictive, see |Wes93] . Under certain assumptions, 
|BMW07l Theorem 2.2] asserts that there exist numerical constants rj,9 G R a, r > and 
^ 1 7 G [—1, 1 + 2-y/or] such that for all s < t < u, 
Q\ ' 

O; /i i\ \t rviTT i (u-t){t-s) ( uX s -sX u X u -X s 
O ; (1.3) Vw[X t \F S:U \ = —— — 1 + r] h 9- 



u(l + as) + r — 7s \ u — s u — s 

(uX s -sX u ) 2 (X u -X s ) 2 (X u -X s )(uX s -sX u ) 
+a — 7 \5 — + T —< \5 ( X - 7)- 



^ | (u — s) 2 (u — s) 2 (u — s) 2 

We will say that a square-integrable stochastic process (X t ) t( z T is a quadratic harness on T with 
parameters (rj, 9, a, r, 7), if it satisfies fll.ip . (11. 2p and (II. 3p on an open interval T C (0, 00). 

Our goal is to construct examples of Markov quadratic harnesses with 7=1 — 2y / o r r. In 
[BMW07, Proposition 4.4], these were called "classical quadratic harnesses. The construction 
follows [BWlOl Section 2] who construct quadratic harnesses with 7 < 1 — l^far from the 
Askey- Wilson integral. Here we use instead some of the generalized Beta integrals from [Ask89] . 

The paper is organized into sections based on the number of parameters in the generalized beta 
integrals. In particular, in Section H] we exhibit explicit transition probabilities for the bridges 
of the hyperbolic secant process, and for completeness in Section [5] we re-analyze the Dirichlet 
process. 
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1.2. Conversion to the standard form. In this section we recall a procedure that we use to 
transform Markov processes with linear regressions and quadratic conditional variances into the 
quadratic harnesses. The following is a specification of [BW09, Theorem 1.1] that fits our needs. 

Proposition 1.1. Suppose (Y t ) is a (real-valued) Markov process on an open interval T C R 
such that 

(1) E(Yf) = a + (3t for some real a, (3. 

(2) For s <t in T, Cov(Y s , Y t ) = M 2 (ip + s){5 + et), where M 2 (tp + t)(5 + et) > on the 
entire interval T , and that 5 — sip > 0. 

(3) For s <t <u, 

(1.4) Var(F 4 |F S) Y u ) = F t>s>u ( X o + Vo ~ - — + ej u ~ Y ' ' {Yu " Ys < 



v It — S M — S (u — s) 2 

where F ttS , u is non-random and xo, #o, ?7o € R ore stzc/i i/ioi X := Xo + c^o + /3#o + /3 2 > 0. 

Denote F t = Y t - E(Y t ). Then there are two affine functions £(t) — M(s—eip) an( ^ m(t) — m(5— e^/i) 

and an open interval T' C (0, oo) such that X t := m(t)Y^t)/m{t) defines a process (X t ) on T' such 
that (11. ip holds and (II. 3p holds with parameters 

(1.5) V = M(6 m + e(2f3 + 9 ))/x, 

(1.6) = M(2/3 + ^ + # )/x, 

(1.7) a = MV/x, 

(1.8) r = M 2 /x, 

(1.9) 7 = l + 2e v / ^- 

Proof. This is |BW09l Theorem 1.1] specialized to x — Xo, V — Vo, — $o, & — 0, r = 1, p = 0, 

a = M, 6 = MV>, c = Me, d = M5. □ 

Remark 1.1. We u>i/Z appZy £/ws only to e — 0, ±1, and xo, #0, £ {0, 1}. 

Remark 1.2. For e < 0, we see that 7 < 1 and r/y/r + 9y/a = M 2 (5 — eip)r}o/x 2 has the same 
sign as r) . 

Remark 1.3. The time domain T' is the image of T under the Mobius transformation t 1— > 

(t + l[>)/(Et + 8). 

Two related transformations are sometimes useful to keep in mind, as they take care of some 
additional non- uniqueness in the final form of (II. 3p . Firstly, if (Xt) is a quadratic harness with pa- 
rameters (77, 9, a, t, 7) then (aXt/ a 2 ) is a quadratic harness with parameters (rj/a, a9, a/ a 2 , a 2 r, 7). 
In particular, if a = and r > 0, then without loss of generality we may take r = 1. And if 
cr, r > then without loss of generality we may take a = r. (So our constructions will lead to 
these two cases only.) 

Secondly, time inversion (tX\/ t ) converts a quadratic harness with parameters (rj, 9, a, r, 7) into 
a quadratic harness with parameters (9, n, r, a, 7), i.e. it swaps the entries within pairs (77, 9) and 
(cr, r). In particular, time inversion maps a quadratic harness with a = 0, r = 1 into a quadratic 
harness with a = 1, r = 0. Similarly, it maps a quadratic harness with parameters a = r and 
rj 2 < 4cr, 9 2 > 4cr into a quadratic harness with parameters a = r and f] 2 > 4er, 9 < 4a. 

2. Four-parameter beta integral 



This section contains the construction of Markov processes based on the four-parameter beta 
integral [Ask89, (8.i)]. After a transformation, these processes become quadratic harnesses with 
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arbitrary a — r G (0, 1), 7 = 1 — 2y/or, and with rj, 6 such that y/rrj + v^"# 7^ 0; parameters 
i], 9 will be required to satisfy also some additional restrictions, of which 1]6 > suffices for 
all constructions to go through. Since the main steps will be repeated several times, first with 
three parameters to cover the case a = 0, and then with two parameters to cover the case 
y/rrj + y/aO = 0, we give here more details so that we can suppress them in the subsequent 
iterations. 

The construction starts with four complex numbers 01,02,03,04 with strictly positive real 
parts. The generalized beta integral |dB72| IWil80] after changing the variable to y/x is 

^ 111=1 ( r (Qj + iy/x)T(aj - Vg) ^ _ 4?r IL<fc<j<4 r K + a i) 



y/x\T{2iy/x)\ T(oi + o 2 + o 3 + o 4 ) 

(2.2) K{a,b,c,d) 



Denote 

T(a + b + c + d) 



4vrr(o + b)T(a + c)T(b + c)T(a + d)T(b + d)T(c + d) ' 

If a, 6, c, d are positive real numbers, or come as one or two conjugate pairs with positive real 
parts, identity ( 12. ip implies that the following function of x > becomes a four-parameter 
probability density function on (0, 00): 

\T(a + iy / x)T(6 + iy/x)T(c + iy/x)T(d + i\/x)\ 2 



(2.3) f(x;a,b } c } d) = K(a } b } c } d)- 



x|r(2i^)| 2 



Proposition 2.1. If a random variable X has density f(x;a,b,c,d) then 

abc + abd + acd + bed 



(2.4) E(X) 
and 

(2.5) Var(X) 



a + b + c + d 

and 

(a + b)(a + c)(b + c)(a + d)(b + d)(c + d) 



(a + 6 + c + d) 2 (a + b + c + d+ 1) 

Proof. The formulas can be read out from the first two orthogonal polynomials [KS98, (1.1.4)], 
but they also follow easily from the formulas 

K(a, 6, c, d) 



E(a 2 + X) =E ((a + iVx)(a - iVx)j 



K(a + 1, 6, c, d) 



and 

a 2 6 2 + (a 2 + 6 2 )E(X) + E(X 2 } = jg ( (fl 2 + x)(fe 2 + x) ) _ ^ [a. b. r. d) 



K(a+ 1,6+ l,c, d)' 

Now using (12.2p and sT(s) = T(s + 1), we get (12.4[) . and after a calculation we get ( 12. 5ft . □ 

Next, we prove a "convolution formula" which will be used to verify the Chapman-Kolmogorov 
equations. 

Proposition 2.2. If m > £/ien 

;>oo 

(2.6) f(y;a } b J c + m } d + m)= / / (y;a,b,m + i^/x,m — iy/x) f(x;a + m,b + m,c,d) dx. 

Jo 

Proof. Re-arranging the factors in ( 12. 3p . we have 

,„m f(x,a + m,b + m,c,d)f (y,a,b,m + iy/x,m - iy/x) . 

(2.7) ■ — r = / (x, m + iy/y, m - iy/y, c, d) . 

j [y, a,b,c + m,d + m) 

Formula (12. 6p now follows, as J °° f (x; fi + iy/y, \i — iy/y, c, d) dx = 1. □ 
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We remark that ( 12 .7p is an analog of |Jam74t (b3)] and will serve similar purposes. Related 
formulas will appear again as (I3.3p . (14. 5p . (I4.14p . and (I5.3p . 

2.1. The auxiliary Markov process. We now define a family of Markov processes (Y t ) t£ T, 
parameterized by A, B, C, D that are either all real and positive or come as one or two complex 
conjugate pairs A = B or C = D, with positive real parts. Without loss of generality we may 
assume that R(A) < 3fc(B) and 9fc(C) < ftp). 

As the time domain for Markov process (Y t ) we take the open interval T = (— 3?(C), §l(A)), 
and as the state space we take (0, oo). We define the univariate distribution of Y t by the density 

(2.8) f t (x) = f(x; A — t, B — t,C + t,D + t). 

For s < t, we define the transition probability C{Y t \Y s = x) by the density 

(2.9) fsM x ) = f(y,A-t,B-t,t-s + iVx,t-s- i^/x) . 
It remains to verify that the above definitions are consistent. 

Proposition 2.3. Formulas (12. 8p and (12. 9p determine a Markov process (Y t )t^T- Furthermore, 
E(Y ) = {ABC + ABD + ACD + BCD) /{A + B + C + D) by ([221) and 

A JD C V) 

(2-10) E(Y t ) = E(Y ) + 2 ATBTcTD t ~ ^ 
For s < t in T , 

(2.11) Cov(r s , Y t ) = M 2 (C + D + 2s)(A + B-2t), 
where 

(2.12) M 2 = (A + C)(B + C)(A + D)(B + D) 

y ' {A + B + C + D) 2 {A + B + C + D + 1) 

In view of (I2.10p . the conditional moments simplify when we express them in terms of 

(2.13) Y t = Y t/2 + t 2 /A, -25R(C) < t < 2M{A), 

with linear mean K(Y t ) — a + fit and the covariance Cov(Y s , Y t ) = M 2 {C + D + s)(A + B — t) 
for s < t. The one-sided conditional moments s <t are: 

(2 ,4) E(?1 | ?lH <4±£^) ? , + ^(<-s) 



A + B-s A + B-s' 
(A + B — t)(t — s) i^A 2 -sA + Y s ^j (B 2 -sB + Y s 

■t\ J -s) — 

For s < t < u in T, 



(2.15) Var(F t |r s ) 

v ' v 1 ' (A + B — s) 2 (A + B — s + 1 



(2.16) E(Y t \Y.X)= { — t)Y - + {t - S)Y " 



u — s 



(2,7, V*m, Y,) = ^MzA 0^ + « 

u — s + 1 \ (u — s) 2 u — s 
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Proof. To verify Chapman-Kolmogorov equations we first use (12. 6p with m = t — s, a = A — t, 
b = B-t,c = C + t, d = D + t. This gives 

poo 

(2.18) f t (y)= / f att (y\x)f 8 {x)dx. 

Jo 

Next we use (I2.6P with m = u — t, a = A — u, b = B — u, c = t — s + iy/x, d = t — s — i\fx to 
verify the Chapman-Kolmogorov equations for the transition probabilities, 

poo 

(2.19) f S)U (z\x) = f S)t (y\x)f t;U (z\y)dy. 

Jo 

Formula (12 .7j) can be now reinterpreted as the formula for the conditional distribution C(Y t \Y s = 
x,Y u = z), given by the density 

(2.20) g(y\x, z) = = f (y- u _ f _|_ u — t — i\fz, t - s + i\/x, t - s - i\fx) . 

Since this is again expressed in terms of the same density (12. 3p . the formulas for the conditional 
mean and conditional variance are recalculated from Proposition 12. II Finally, we use (I2.10p . and 

(2.21) Var(y t ) = M 2 (A + B - 2t)(C + D + 2t), 

which is calculated from (12. 5p . and (I2.14p . to compute K(Y s Y t ) and we get (12. lip . □ 

Corollary 2.4. (Y t ) can be transformed into a quadratic harness with covariance (II. ip and the 
conditional variance (II. 3p with parameters 

(A + B + C + D) 2 



(2.22) 


7/ 


+ 9 


(2.23) 





-7] 


(2.24) 


a 


= T 


and 7 = 1 


- 2^ 


far 



y/(A + C)(B + C)(A + D){B + D){A + B + C + D + 1) ' 

(C - D) 2 - (A - Bf 

yJ(A + C)(B + C)(A + D)(B + D)(A + B + C + D + 1) ' 
1 



Proof. We apply Proposition 11.11 with parameters 

" = A + B + C + D ^= A + B + C + D > e = - 1 ^ = C + D ' 6 = A + K 

The only non-zero parameters in (11.4p are 7]q = tq = 1. □ 
Remark 2.1. TTie quadratic harness is defined on the interval 



= f C + D-2^(C) C + D + 2$(A) \ 
= \ A + B + 25R(C) ' A + B - 2^{A) J ' 



In particular, T' = (0, oo) if A = B and C — D. It is plausible that by allowing transition 
probabilities and univariate laws with discrete components, this interval could be extended to 
(0, oo) in all cases when dt(A + B) > and dt(C + D) > 0. 
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2.2. The admissible range of 77, 9. In this section we study which collections of parameters 
correspond to quadratic harnesses from the previous construction. Given a = r, 7 = 1 — 2a and 
77, 9 such that i] + 9 ^ 0, without loss of generality we may assume that 77 + 9 > 0. For if we can 
find a quadratic harness (X t ) for one such set of parameters, then (—X t ) is a quadratic harness 
on the same time domain, with the same a, r, 7, but with —77, — 9 instead of 77, 9. 

Once we restrict ourselves to the case 77+$ > 0, we want to know for which 77, 9, a = r, 7 = 1— 2a 
we can find A,B,C,D that satisfy the equations from Corollary 12 .41 and satisfy the constraints for 
the construction of the Markov process (Y t ). We will see that we can always find such A, B, C, D 
if either i]9 > (which under the assumption 77 + 9 > 0, is equivalent to 77 > 0, 9 > or vice 
versa) or when the sign of t]9 is arbitrary but ry 2 < 4a and 9 2 < 4t. 

To proceed, we rewrite the equations from Corollary 12.41 in equivalent form: 

(2.25) A + B + C + D = (I -a) /a, 

(2.26) (A + C)(B + C)(A + D)(B + D) (i ~' Ti 

(2.27) (C - D) 2 - (A - B) 2 



(77 + 9) 2 a 3 ' 
(9-y)(l-a: 
(77 + 9)a 2 



2 



2.2.1. Hyperbolic case. We first show that quadratic harnesses exist for any 77, 9 such that r] + 9 > 
0, rf < 4a and 9 2 < 4a. This is because in this case the system of equations fl2.25ti2.27j) is solved 
by two conjugate pairs A = B and C = D with A, C given by 

A = + 1(1 C = ^- ®(A) + * (1 7 )v/ ^ 

V ; 2(77 + 0)0- 2a V ; 2{rj + 9)a 

with arbitrary < < 

The apparent non-uniqueness in this solution and in others is in fact illusory, as it corresponds 
to the translation of the time domain T. This translation does not affect neither the transition 
probabilities of the final quadratic harness, nor the final time domain, which by Remark 12.11 is 
T' = (0,oo). 

2.2.2. Next we go over the remaining choices of pairs {f],9), and confirm that in each case we 
can always find a quadratic harness when t]9 > 0. 

We first consider 77, 9 such that 77 + 9 > 0, r\ 2 < 4a and 9 2 > 4a. Then quadratic harnesses 
exist if 4a + 77 2 + 2776* > 0. Indeed, in this case the system of equations fl2.25H2.27j) is solved with 
one conjugate pair A = B. The solutions are 

2(77 + 9)a 

(77 + 9 + V9 2 - 4a) (1 - a) 
D = ^ — -L± - U(A) . 

The restriction 4a + rj 2 + 2rj9 > guarantees that 9 2 — 4a < (9 + t?) 2 so one can find dt(A) > 
such that C > 0; then D > follows automatically. 

The restriction 4a + rj 2 + 2i]9 > holds, in particular, if 77 > 0, as then 4a + rj 2 + 2?7# > 
2r] 2 + 2r]9 = 277(77 + 9) > 0. 
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We remark that the left endpoint of the time domain T' here is 0, see Remark 12.11 This is 
of interest, since for such domains the one-sided conditional moments (I2.14p and (12.151) imply 
uniqueness of the quadratic harness. 

Finally, if rj + 9 > are such that rj 2 > 4er and 9 2 > 4cr, then under the condition rj + 9 > 
yV -4a + V9 2 - 4<r one can choose a small enough A > such that 



(77 + 6)a 

(rj + 9- J if - 4a - y/9 2 - 4a) (1 - a) 
2( V + 9)a 

(rj + 9- Jrf - 4a + V9 2 - 4a) (1 - a) 

D = 7 ^ A 

2(r] + 9)a 



are all positive. In particular, if r], 9 > then r\ > \Jri 2 — 4a and 9 > y/9 2 — 4a so the above 
solution will indeed give us a quadratic harness on a finite interval T' . 

3. Three parameter beta integral 

For a > and b, c real positive or a complex conjugate pair with positive real part, define the 
following density on [0, oo): (See |Ask891 (7.i)] or jKSMl Section 1.3]) 

|r(a + iy/x)T(b + iy/x)T(c + iy/x)\ 2 
1 j 9[X] a ' ' C) ~ 4nT(a + b)T(a + c)T(b + c)v^|r(2z^)| 2 ' 

As previously, it is straightforward to use properties of the gamma function to get formulas for 
the mean \x and the variance a 2 , 

(3.2) /i = ab + ac + be, a 2 = (a + b) (a + c)(b + c). 
The relevant version of (12.71) is 

io ?\ g(x;a + m,b,c)g(y;a,m + iy/x,m-iy/x) . 

(3.3) — ■ r = f{x,m + i^y,m-t^y,b,c), 

g(y; a,o + m,c + m) 

where x,y,m > 0. 

Let A G R and let B, C be either real or a complex conjugate pair, and without loss of 
generality we assume that in the real case B > C . Suppose in addition that A + dt(C) > 
so that T = (— 9ft(C),v4) is non-empty. Then from (13.31) we get again a Markov process (Y t ) teT 
with univariate distributions on the state space (0, oo) defined by the densities g(x; A — t, B + 
t,C + t), with transition probabilities defined for s < t in T and x,y > by the densities 
g(y; A — t,t — s — i^/x,t — s + i^/x), and whose two-sided conditional laws are again given by 
Wilson's density (I2.20p . In particular, after we make substitution (I2.13P formulas (I2.16P and 
(I2.17P for the two-sided conditional mean and variance hold. 

As previously, parameters A, B, C affect only the mean and the covariance of (Yt): 

(3.4) E(K t ) = -t 2 + 2 At + AB + AC + BC, V&r(Y t ) = (A + B)(A + C)(B + C + 2t) . 
Passing to the centered process (12.131) . the one-sided conditional moments are: 

E(Y t \Y s ) = A(t -s) + Y s , Var(F t |y s ) = (t - s) (a 2 -sA + Y s ' 

In particular, the above formula for E(Y^|y s ) gives 

Cov(r;, Y t ) = (A + B)(A + C)(B + C + mm{t, s}) . 
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Then the transformation from Proposition 11.11 takes a particularly simple form. Markov process 

Y t _ B _c - E(y t _ B -c) 



y/(A + B)(A + C) 



X t = 

is a quadratic harness with parameters 

(3.5) r] = 1 , 

y/{A + B){A + C)' 

3.6 = 

(3.7) <7 = , 

(3.8) r = 1. 

This gives us a family of quadratic harnesses arbitrary positive values for parameters 77, 9, with 
r = 1, a — 0. Other values of parameters are now produced by routine transformations that 
were mentioned in the introduction. To swap the roles of a, r one uses time inversion (tXut). 
Taking (— X t ) we get arbitrary negative values of 77, 8, covering all possible non-zero values of 
the same sign (t]6 > 0). Finally, transformation (X at / \fa) produces arbitrary positive values for 
parameter r. 

Remark 3.1. The quadratic harness is defined on 

T' = (9fi(C-£),oo). 

In particular, T' = (0, 00) if B = C. It would be interesting to see if the construction could be 
modified to yield T' = (0, 00) also for real B ^ C . 

Remark 3.2. Formula (13. 3p indicates that bridges of the three-parameter quadratic harnesses 
with o = are the (transformations of) four-parameter quadratic harnesses from Corollary \2. 4 
It would be interesting to see if this holds true also in the cases without densities. 



4. Two-parameter beta integral 

According to |Ask89j (5.i)], see also [KS98i Section 1.4], the following is a probability density 
on R when c = a,d = b have positive real part. 

T(a + b + c + d)T (a + ix) T (b + ix) Y (c — ix) Y id — ix) 
(4.1) ip{x\ a, b, c, d) = — 



2vrr(a + c)Y(b + c)Y(a + d)Y(b + d) 
The analog of Proposition 12.11 is 
Proposition 4.1. If a random variable lei has density <p(x; a, b, c, d), then 

&(a)S(6) +3fc(6)9f(a) 



(4.2) E(X) 



(4.3) Var(X) 



tt(a + b) 

9t(a)9t(6) ((9ft(a + b)f + (S(a - b)f) 
(K(a + 6)) 2 (2K(a + 6) + l) 
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Proof. Denote by K(a, b, c, d) = 2 nr(a+c)T(b+c)r(a+d)r(b+d) tne normalizing constant in Then 

POO 

(4.4) / x(p(x;a,b,c,d)dx 

((a + ix)(c — ix) — (b + ix)(d — ix) + bd — ac) (f(x; a, b, c, d) dx 



i(c + b — a — d) 



K(a,b,c,d) ( 1 1 , , \ i(ab — cd) 



i(c + b — a — d) \K(a + l,b,c + l,d) K(a, b + 1, c, d + 1) / a + 5 + c + d 

Substituting a = 9fc(a) + z3(a), & = &(&) + c = »(a) -z3f(o), d = »(&) we get (Ojl . 

The variance comes from a similar calculation: 

E(X 2 ) - (E(X)) 2 = K(a + 1, 6, c + 1, d) - (c - a)E(X) - ac - (E(X)) 2 

(a + c)(& + c)(a + d)(b + d) 



(a + 6 + c + d) 2 (a + & + c + d+l)' 

□ 



The analog of Proposition 12.21 is based on the identity 

(p(y; a, m — ix, a, m + ix)<p(x; a + m, b, a + m, b) - 

(4.5) — ~- = tp{x;b,m-iy,b,m + iy). 

(p{y; a,b + m,a,b + m) 

Thus, given complex parameters A, B, such that 3t(A + B) > 0, let T = (— 3t(B), 31(A)). For 
s < t in T, the univariate densities on R, 

(4.6) f t (x) = <f(x, A — t, B + 1, A — t, B + 1) , 
and the transition probabilities 

(4.7) f St t(y\x) = (p(x, A — t,t — s — ix, A — t,t — s + ix) , 

satisfy the Chapman-Kolmogorov equations ( 12. 18ft and ( 12.191) . Let {Y t )teT denote the corre- 
sponding Markov process. Then from ( 14. 2 ft and ( 14.31) we get 

(4 . 8) EW) = jgz| ( _ »|B , v ar ,y«) = *W) - oot + 1) , 

where 

(3(A- J B)) 2 + (K(A + J B)) 2 



(4.9) M 

y ' (3t(A + B)) 2 {23t{A + B) + l) 

Since (14. 2p also gives 

for s < t, from (14.81) we further calculate 

(4.10) Cov(r s , Y t ) = M 2 (^(A) - t)(3t(B) + s). 

Next we compute conditional moments. For s < t < u in, the two-sided conditional density of 
C(Y t \Y s = x, Y u = z) is given by 

g(y\x, z) = <p(y; t — s — ix,u — t — iz,t — s + ix,u — t + iz) . 

So from (jOl . 

E(y.|r.,y.) = (- < ) y - + (t -) y v 
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and from (14. 3[) we get 

(t-s)(u-t) / , (Y u — Y s ] 2 



(4.11) Var{Y t \Y s ,Y u 

From Proposition O applied with Xo = 1, a = - ^g'^ . f> = 



(2m - 2s + 1) \ (u - s) 2 



2K(A + 5) + 1 
2(9f(A-S)) 



vPP + 5) + 1) pp - 5)) 2 + (ip + 5)) 2 ) 

1-CT 



From the first equation, we see that $t(A + B) = The second equation determines ^(A — £?) 
as a real number iff 9 2 < At. This proves the following. 

Proposition 4.2. For every a G (0, 1) and r\ G (— 2a/cx, lyfo), there is a quadratic harness on 
(0, oo) with parameters rj,9 — —T), a, r = a, 7 = 1 — 2a. 

4.1. Bridges of the hyperbolic secant process. Bridges of all Meixner processes are de- 
scribed in [BW091 Proposition 4.2 and Remark 4.1]. According to these results, bridges of 
Meixner processes are quadratic harnesses with r\\Jr + 6y/a~ = 0. When or > 0, then depending 
on the sign of 9 2 — At, such processes arise as bridges of negative binomial, gamma, or the hy- 
perbolic secant process. In [BW09| the bridges of hyperbolic secant process were not described 
explicitly, so we identify their transition probabilities here. 

The following integral is due to Meixner [Mei34[ page 13], and is listed as [Ask89l (4.i)]: 



(4.12) 



loo' 1 M (2cosf) 2 « 



The integral is well defined for real a > and — n < (3 < ir. Denote by f(x; a, (3) the corresponding 
density, i.e. 

f2cos £) 2a 

(4.13) f(x; a, 0) = ^— jL |r( a + zx)| V* , 

and by X the corresponding random variable. 

Differentiating (I4.13p with respect to (3 and integrating the answer we get E(X) = atan(f) 
and Var(X) = |asec 2 (|). It is known that the corresponding Markov process has independent 
increments: the univariate law of Y t has density ft{x) = f(x; A — t, (3) and the transition densities 
are f s j(y\x) = f(y — x; t — s,(3). One can verify also Chapman- Kolmogorov equations directly 
from the analog of Proposition 12.21 which is based on the identity 

(4.14) ■ = ip(x;a,m-iy,a,m + iy). 

f(y;a + m,/3) 

The right hand side of ( I4.14p integrates to 1 because of (14. ip . (This gives an "elementary" proof 
of the well known fact established by Laha and Lukacs |LL60t Lemma 2] that the hyperbolic 
secant laws form a convolution semi-group.) 

The following proposition describes in more detail bridges mentioned in [BW091 Remark 4.1]. 
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Proposition 4.3. All bridges of a hyperbolic secant process are described by formulas (14. 6p and 
(14. 7p . Conversely, all quadratic harnesses with 0< err <1 ; 7 = 1— and w,9 G R such 

that r]y/r + Q^fo = and 6 2 < At can be realized as such bridges. 

Proof. From (I4.14p we see that for a hyperbolic secant process (Yj), the two-sided conditional 
law of £(Y t \Y s = x, Y u = z) is given by 

(4.15) g(y\x, z) = ip(y — x;t — s,u — t — %{z — x),t — s,u — t + i(z — x)). 

Inspecting formula (14. ip . we see that 

<p(y — x;t — s,u — t — i(z — x),t — s,u — t + i(z — x)) =(p(y;t — s — ix,u — t — iz,t — s + ix,u — t + iz) 

= (p(y; u — t — iz,t — s — ix,u — t + iz,t — s + ix) . 

So identifying this with the univariate law of the bridge at time S < t < U, conditioned at S < U, 
we can read out that the bridge corresponds to the Markov process with transition probabilities 
flUD , where A = U - iY v , B = -S - iY s . 

□ 

5. Standard beta integral 
In this section we use the well known beta density 

(5.1) f(x ] a } b) = ^±^ ) x a - 1 (l-x) b - 1 

to re-derive the quadratic harness properties of the one-parameter family of Dirichlet processes 
from [BW09, Eample 4.1]. (Here, the density is on < x < 1, and the parameters satisfy 
a,b > 0.) It is well know that the corresponding random variable X has moments 

(5.2) E{X) = -?—, and Var(X) = . 

v ' v ' a + b v ' (a + b) 2 (a + b + l) 

The analog of (12.71) is the algebraic identity 

r^— f (t-^-'i Tn, b) f (x; a, b + m) \ 

(5.3) l - xJ 1 M ' = -/ s/y; «, m). 

f{y;a + m,b) y 

In particular, we have a "convolution formula", 

f y 1 fy-x \ 

(5.4) / / I ; m, 6 I f(x; a,b + m)dx = f(y; a + m,b). 







1 — x V 1 — x 



Given A > 0, we now use (15 .4p to define the Markov process (Y t )o <t< A by specifying its univariate 
laws as 

ft(x)=f(x]t,A-t), 16 [0,1], 
and for s < t, y > x its transition probabilities as 

f.Mx) = ^—f (^;t-s,A-t 



1 — x \1 — x 

A calculation based on (15 .4p shows that these expressions indeed satisfy the Chapman- Kolmogorov 
equations, so Markov process (1^)^(0,^) is well defined. (The same conclusion can be reached via 
probabilistic arguments, as these processes arise as bridges of the gamma process.) 
From (JE2D, E(Y t ) = t/A and Var(F t ) = ^a+i) » and with 

some more work one can read out 

that Cov(F s , Y t ) = -^^y for s < t. 
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Since f 15 . 3 j) shows that two-sided conditional laws are also beta, from (15. 2 p we can read out the 
conditional moments 

E(Y t - Y S \Y S , Y u ) = —(Y u - Y s ), 
u — s 

var(y t |y s ,y u ) = {t ~^ u ~'\ u (Y u - Y s f. 

[u — s) z ((u — s) + 1) 

Applying Proposition 11.11 with M = A J^+i > P = V^j = A, e —1 (the remaining parameters 
are 0), we see that (Y t ) can be transformed into a quadratic harness on X" = (0, oo) with 
parameters 

rj = -9= -2/vTTA a = t = 1/(1 + A), 7 = 1- 2/(1 + A) 2 . 
(This is consistent with [BW09, Eample 4.1].) 

Acknowledgement. We would like to thank Arthur Krener and Ofer Zeitouni for information 
on reciprocal processes, and to J. Wesolowski for several related discussions. This research was 
partially supported by NSF grant #DMS-0904720, and by the Taft Research Center. 

References 

[Ask89] R. Askey. Beta integrals and the associated orthogonal polynomials. Number Theory (ed. K. Alladi). 

Lecture Notes in Mathematics, 1395:84-121, 1989. 
[BMW07] Wlodzimierz Bryc, Wojciech Matysiak, and Jacek Wesolowski. Quadratic harnesses, q-commutations, 

and orthogonal martingale polynomials. Trans. Amer. Math. Soc, 359:5449-5483, 2007. 

arxiv.org/abs/math.PR/0504194. 
[BW09] Wlodek Bryc and Jacek Wesolowski. Conditioning of quadratic harnesses. 

http://arxiv.org/abs/0903.0150, 2009. 
[BW10] Wlodek Bryc and Jacek Wesolowski. Askey- Wilson polynomials, quadratic harnesses and martingales. 

Annals of Probability, 38:1221-1262, 2010. arxiv.org/abs/0812.0657. 
[dB72] L. dc Branges. Tensor product spaces. J. Math. Anal. Appl, 38:109-148, 1972. 

[Jam74] Benton Jamison. Reciprocal processes. Z. Wahrscheinlichkeitstheorie und Verw. Gebiete, 30:65-86, 
1974. 

[JP88] Jean Jacod and Philip Protter. Time reversal on Levy processes. Ann. Probab., 16(2):620-641, 1988. 

[KS98] Roelof Koekoek and Rene F. Swarttouw. The Askey scheme of hypergeometric orthogo- 
nal polynomials and its q-analogue, 1998. Delft University of Technology Report no. 98-17, 
http://fa.its.tudelft.nl/ koekoek/askey.html. 

[LL60] R. G. Laha and E. Lukacs. On a problem connected with quadratic regression. Biometrika, 47(300) :335- 
343, 1960. 

[Mei34] J. Meixner. Orthogonale polynomsysteme mit einer besonderen gestalt der erzeugenden funktion. J. 

London Math. Soc, 9(6), 1934. 
[MY05] Roger Mansuy and Marc Yor. Harnesses, Levy bridges and Monsieur Jourdain. Stochastic Process. 

Appl, 115(2):329-338, 2005. 
[Wes93] Jacek Wesolowski. Stochastic processes with linear conditional expectation and quadratic conditional 

variance. Probab. Math. Statist., 14:33-44, 1993. 
[Wil80] J. A. Wilson. Some hypergeometric orthogonal polynomials. SIAM Journal on Mathematical Analysis, 

11:690, 1980. 



Department of Mathematical Sciences, University of Cincinnati, PO Box 210025, Cincinnati, 
OH 45221-0025, USA 



